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Abstract: We carry out an analysis of the effect of a quenched magnetic field on 
the dynamics of the spherical Sherrington-Kirkpatrick spin-glass model. We show 
that there is a characteristic time introduced by the presence of the field. Firstly, 
for times sufficiently small the dynamic scenario of the zero field case - with 
aging effects - is reproduced. Secondly, for times larger than the characteristic 
time one sees equilibrium dynamics. This dynamical behaviour is reconciled with 
the geometry of the energy landscape of the model. We compare this behaviour 
with experimental observations at a finite applied field. 
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1 Introduction 



For some time there has been much interest in aging phenomena and in general 
non-equilibrium behaviour jjj]. For a wide range of models one finds that if the 
thermodynamic limit is taken before the long-time limit, for quite generic initial 
configurations of the system, one never sees an equilibrium regime characterised 
by time translational invariance of the correlation function and the satisfaction 
of the fluctuation dissipation theorem [0, ||. Recently the authors presented a 
detailed examination of the dynamics of the spherical Sherrington-Kirkpatrick 
(SK) spin-glass model || in Ref.||, hereafter refered as I. In general the initial 
conditions do not lead the system to equilibrium and the correlation function 
was shown to exhibit the aging phenomenon. This model joins a whole range of 
models where non-equilibrium behaviour is present, however the precise physical 
mechanisms responsible for aging in various models can be completely different 
but still give the same mathematical behaviour for various dynamical functions. 
Aging can arise due to the existence of large energy barriers in the system (see e.g. 
Refs.[[|, |5]), also due to the existence of zero modes in the free-energy landscape 
as in the case of the models considered in Refs.[^|, 0], due to a combined effect 
of large energy barriers and extended flat regions in the energy landscape as in 
mean-field models 0, U and also via domain growth or coarsening mechanisms in 
the system ||. In all these cases the correlation functions exhibit aging effects. 
However, the domain growth as well as the pure zero-mode mechanisms do not 
seem to be enough to capture aging in quantities such as the magnetisation. The 
effects of small perturbations in such models are rapidly erased and henece the 
response decays too fast to reproduce the slow relaxations of glassy systems. 

In this paper we revisit the spherical SK model but with the addition of a 
quenched magnetic field. For simplicity we work with a random magnetic field 
but this does not change the underlying physics of the problem. The theory of 
linear response applied to this model would suggest that, for sufficiently small 
fields, the behaviour of the system is not drastically altered and that one should 
recover the type of aging phenomena observed in I. However the random field 
strength appears coupled to a nonlinear term in the equation for the dynamical 
Lagrange multiplier imposing the spherical constraint and hence it is not apparent 
that linear response theory should hold. We shall concentate on the case of zero 
temperature where the discussion of both the statics and dynamics is most clear. 
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This is however sufficient to bring out the sentient points of the problem. One 
finds that the random magnetic field introduces a time scale into the problem. 
Below this characteristic time one indeed sees aging phenomena and the effects 
seen in I are reproduced. Above this timescale the system does indeed reach 
an equilibrium state. With the benefit of this insight we carry out a simple 
analysis of the energy landscape and point out that it is the disappearance of the 
preponderance of zero modes in the energy landscape which is responsible for the 
ultimate equilibrium behaviour, adding weight to the assertions made in I. 



2 Dynamics 

We shall recall briefly the definition of the model. The Hamiltonian is given by 

H = --J2JijSiSj-J2 h i s i> (!) 

ij i 

subject to the spherical constraint s l — (N being the number of spins). 
The matrix J is a random symmetric matrix with independently Gaussianly dis- 
tributed components scaled with N to give the Wigner semi-circle law distribution 
for the eigenvalues \x in the thermodynamic limit, that is 

p( yU ) = i- v /4^, 6 [-2,2]. (2) 

The terms hi represent the random field which is chosen to be Gaussian with zero 
mean and variance h 2 . The dynamical equation for the evolution of the spins is 
then given by (at zero temperature) []3], |10| 

Os ■ 

= JijSj - z(t)si + hi. (3) 

We shall assume uniform initial conditions for the spins and proceed by diag- 
onalizing the equations of motion and computing the dynamical Lagrange mul- 
tiplier z self-consistently. Defining 

Q(t) = exp (7* z{t')dt^j (4) 
we find it obeys the following equation 

n 2 (t) = fit) + h 2 j* dt' J* dt"f (t - t -Af- ] j WMt"), (5) 
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where 

f(t) = J fMdpP*. (6) 

In constrast with the zero field case the equation satisfied by the dynamical 
Lagrange multiplier is now non-linear (where as before it was a linear Volterra 
equation 0). The equation @ is evidently very difficult to solve, however it is 
clearly causal and accessible to numerical solution. We precede by making the 
following ansatz on the asymptotic form for Q 

fi(t)~ce At , (7) 

where c is some positive constant. Futhermore if we assume a priori that A > 2 
then because f(t) ~ e At (2t)^ / y/Air for sufficiently large times we may make the 
asymptotic approximation 

n 2 (t) = h 2 f df f dt"f ( t - t -^-] n(t')n(t") . (8) 



Assuming the Q is well behaved and bounded for small times we may also as- 
sume that the exponetial behaviour dominates the double integral to obtain the 
following equation determining A: 

The equation (|D yields the solution 

2 + h 2 

A = . 10 

VTTh 2 

hence we see a posteriori that indeed A > 2 for h > 0. Of course the fact that we 
have found an asymptotic solution for Eq.(§) does not mean that it is the solution 
that matches with the given initial conditions (the problem of matching in these 
kinds of systems is in general outstanding, e.g. see [0). However we have checked 
numerically that this is indeed the case. The results of a numerical integration 
of Eq.(^|) for h = 0.5 and h = 0.4 are shown in figures 1 and 2 respectively. The 
function plotted is actually u(t) = e~ 2t Q(t). The figures clearly show the onset of 
an exponential behaviour of u(t) around the characteristic time scale r(h), with 
the analytically predicted rate (plotted as the dashed line) ; also shown is the zero 
field result (dotted line). 
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Hence we see via the above analysis that the addition of the random magnetic 
field introduces a time scale r into the problem such that for times sufficiently 
large compared to r the contribution of the first term of the righthandside of 
Eq.(|5|) becomes negligible compared to the second term. (It is in fact the first 
term that contains all the information about the initial condition; the fact that it 
becomes negligible demonstrates that the initial condition is completely forgotten 
at large enough times.) The analysis in addition shows that 

^(A-2H(^-2)~\ (11, 

which clearly diverges when h — > as h~ 2 . However for sufficiently small fields 
and/or times it is the first term of the righthandside of Eq.(^]) which dominates, 
and this leads to the aging phenomena exhibited in || within these field/time 
regimes. 

The explicit demonstration of equilibrium behaviour for sufficiently large 
times is now rather trivial. The correlation function is given by 



/ ( t±f] + h 2 fds f ds'f ( t + t ' s s ' ) n(s'Ms) 



(12) 



Inserting the asymptotic form for Q simply yields 

C(M')~1, (13) 

for sufficiently large times. Hence the system ultimately reaches its equilibrium 
state, which is simply the point of minimal energy in the zero temperature case. 
The energy density £(t) can be shown to obey 

S(t) = -~]og(a)-^a(t), (14) 
where a is the average spin/field correlation 

a(t)= ^ EM(W= _|_/; dt7 (^ fi(0 . (15) 

One finds that 

£{t) -> -VTTh 2 (16) 



and 



a(t) -> ; = (17) 
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3 Zero Temperature Statics 



In this section we carry out a simple analysis of the geometry of the energy 
landscape for the model in the presence of a random magnetic field. The energy 
density corresponding to a spin configuration s is given by 

NS = --s • Js - h ■ s + - (s • s - N) , (18) 

where A is the (static) Lagrange multiplier imposing the spherical constraint. The 
variational equations yield the stationary solutions 

s = (A-J)- 1 h (19) 

where A satisfy the equation 

ihtwrw*/^^. (20) 

For A e [—2, 2] the integral above diverges, and hence for finite h we must look 
for solutions outside this interval. There are only two and are given by 

The corresponding values of the energy density can be shown to be 

^-¥'JA»-^ < 22 > 

The minimum of the energy is given by A + , and A_ gives the maximum energy 
(which must also exist according to Morse theory). The minimum energy and 
the static value of a are in agreement with the dynamically calculated limits in 
equations fll6|) and (|17D. 

If one analyses the zero magnetic field case one finds that there are N sta- 
tionary points for the energy, each one corresponding to an eigenvector of the 
interaction matrix. Hence the effect of the field statically, has been to reduce the 
number of stationary points from a macroscopic number down to just two. (This 
is in sharp contrast with the reduction of the number of metastable states in the 
SK model by a magnetic field, where the number always remains macroscopic for 
finite fields ||11|| .) Therefore from almost every starting configuration (except the 
maximum which is however measure zero), the system has a drift towards the 
equilibrium state. It was postulated in || that the large number of flat regions in 
the zero field landscape was responsible for the aging phenomena observed. Here 
we see clearly that once these are erradicated an equilibrium state is achieved. 
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4 Conclusions 



We have analysed the influence of a magnetic field on the dynamics of the spher- 
ical SK spin-glass. At variance with the prediction of linear response, the aging 
behiour of the system is completey erradicated after a characteristic time which 
decreases with increasing field strength. The failure of linear response is presum- 
ably due to the inherent nonlinearity introduced in the equation of motion for 
the dynamical Lagrange multiplier. However for sufficiently short times, linear 
response is presumably valid and one does indeed see that below this time scale 
the aging phenomenon is present. In fact this model explicitly demonstrates the 
phenomenon of interrupted aging discussed in [|12| . 

Also we have related the ultimate arrival at equilibrium of the system to the 
geometry of the energy landscape of the model, showing that the field reduces the 
number of stationary (flat) points in the energy landscape from a macroscopic 
number to just two. This is in agreement with the notion that in the spherical SK 
spin-glass in the absence of a random magnetic field, it is the zero modes in the 
system which give rise to the aging phenomenon rather than tunnelling through 
large energy barriers. 

The behaviour of the spherical SK model is reminiscent of what is predicted 
by the scaling approach to spin-glasses |14|]. Within the droplet picture, the low- 
temperature phase is characterised by only two pure states as for the spherical 
SK model. The other hall-mark of the droplet picture is the absence Jl3| of a 
genuine h — T transition line (AT line) |15|] and hence the destruction of the 
thermodynamic spin-glass phase by any applied field. However, a pseudo AT line 
defined as the locus of points at which the system gets frozen on experimental 
times scales moving down with increasing time scales is claimed to exist in this 
approach For the spherical SK model, statically there is no AT line but 

dynamically we find that the non-equilibrium aging behaviour exists below a 
transition line that moves towards zero field when the times explored increase. 

Experimentally there is no clear evidence for the appearance of an equilibrium 



regime for spin-glasses |16| under the influence of a magnetic field. However, 
we expect that the behaviour observed in this letter may be present in other 
disordered systems such as the problem of a driven particle in a random potential. 
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List of Figure Captions 

Figure:l Graph of u>(t) verses time for h = 0.5. 
Figure:2 Graph of u>(t) verses time for h = 0.4 



Figure: 1 . 




